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Abstract 


The  report  summarizes  a  research  program  focused  on  two  major  issues  related  to 
the  achievement  of  high  Reynolds  number  quiet  flow  in  hypersonic  tunnels,  namely, 
1)  automated  optimal  design  of  quiet  hypersonic  tunnels,  and  2)  tunnel  startup  in 
the  presence  of  blunt  models.  In  the  former  case,  a  fully  automated  optimal  design 
methodology  was  developed  to  determine  the  optimal  shape  of  the  supersonic  nozzle 
to  aehieve  laminar  flow  on  the  nozzle  walls  and  maximize  the  length  of  the  region 
of  quiet  flow  in  the  test  section.  In  the  latter  ease,  detailed  time-accurate  numerical 
simulations  were  performed  for  a  eone  model  of  different  angles  in  the  Purdue  Quiet 
Tunnel.  Results  demonstrated  the  limiting  size  of  the  model  that  permitted  startup  of 
the  tunnel. 
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1  Introduction 


Accurate  prediction  of  the  location  and  extent  of  boundary  layer  transition  in  hypersonic 
flight  is  both  a  critical  parameter  in  hypersonic  vehicle  design  and  an  unsolved  problem.  The 
location  of  transition  affects  estimates  of  aerodynamic  heating,  skin  friction  drag  and  other 
boundary  layer  properties.  Transition  experiments  have  been  performed  in  conventional  wind 
tunnels  for  decades;  however,  most  of  the  experimental  data  obtained  from  these  facilities 
are  contaminated  by  the  high  levels  of  noise  that  radiate  from  the  turbulent  boundary  layers 
normally  present  on  the  nozzle  walls.  These  high  noise  levels  can  cause  transition  to  occur 
at  Reynolds  numbers  (based  on  arc  length)  which  are  an  order  of  magnitude  earlier  than  in 
flight  [4],  and  can  also  affect  parametric  trends  for  transition  [111. 

Experimental  hypersonic  transition  research  in  ground-based  facilities  requires  quiet  flow 
in  the  test  section.  This  implies  laminar  flow  on  the  test  section  walls.  Schneider  [13] 
reviewed  the  history  and  status  of  quiet  flow  hypersonic  wind  tunnels.  Beckwith  [5,  6]  led 
the  design  and  fabrication  of  quiet  tunnels  at  NASA  Langley.  A  quiet  Mach  3.5  tunnel  was 
successfully  demonstrated  at  NASA  Langley  in  the  early  1980s,  followed  by  a  quiet  Mach  6 
hypersonic  facility  in  the  mid-1990s.  The  Boeing/AFOSR  Mach  6  Quiet  Tunnel  (BAM6QT) 
at  Purdue  University  was  constructed  during  1995-2001.  It  is  designed  to  achieve  quiet  flow 
for  stagnation  pressures  up  to  150  psia  (1000  kPa),  corresponding  to  a  unit  Reynolds  number 
of  3.4- 106  ft-1  (11  •  106  m_I)  for  a  stagnation  temperature  of  433  K  at  Mach  fi.  It  is  presently 
the  only  operational  hypersonic  quiet  tunnel  in  the  world  [14]. 

The  conventional  approach  to  design  of  the  supersonic  nozzle  for  quiet  hypersonic  tunnels  is 
trial-and-error.  An  initial  design  is  obtained  using  the  conventional  method  of  characteris¬ 
tics  [16]  based  upon  an  assumed  length  for  the  test  section.  The  displacement  thickness  of  the 
assumed  laminar  boundary  layer  on  the  nozzle  wall  is  subtracted  from  the  inviscid  contour 
to  define  the  actual  shape.  A  linear  stability  or  parabolized  stability  analysis  is  performed 
for  the  laminar  boundary  layer  on  the  nozzle  wall  [7].  The  location  of  transition  is  defined  by 
the  growth  factor  N  in  the  stability  analysis.  The  Mach  cone  emanating  from  this  location 
(assuming  a  circular  cross-section  for  test  section)  and  the  Mach  cone  emanating  from  the 
inviscid  method  of  characteristics  form  a  back-to-front  double  cone  region  which  defines  the 
quiet  flow  region  in  the  test  section.  This  is  illustrated  in  Fig.  1  for  the  BAM6QT.  This 
procedure  is  then  repeated  until  the  desired  size  of  the  quiet  flow  region  in  the  test  section 
is  achieved.  The  obvious  weakness  of  this  approach  is  the  use  of  trial-and-error  to  determine 
the  optimal  quiet  flow  region. 

The  first  objective  of  the  research  program  is  the  development  of  an  automated  optimal 
design  methodology  for  quiet  hypersonic  tunnels.  The  trial-and-error  approach  is  replaced 
with  an  automated  methodology  which  determines  the  maximum  quiet  flow  test  section  size. 
The  methodology  is  described  in  Section  2. 

An  important  concern  for  all  quiet  hypersonic  tunnels  is  starting  the  tunnel  in  the  presence 
of  a  test  model.  Although  slender  vehicles  are  the  primary  concern  in  many  transition 
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Figure  1:  Test  Section  of  Boeing- AFOSR  Mach  6  Quiet  Tunnel 

experiments,  blunt  vehicles  are  also  affected  by  transition  [12].  Shooks  emanating  from  the 
nose  of  the  test  cone  and  bow  shocks  from  blunter  models  interact  with  the  boundary  layer 
on  the  tunnel  wall.  While  disturbances  in  supersonic  flow  can  only  travel  downstream, 
disturbances  in  the  subsonic  boundary  layer  flow  in  the  test  section  can  lead  to  separated 
flow  upstream  in  the  tunnel  nozzle  [17].  A  laminar  boundary  layer  is  more  likely  to  separate 
than  turbulent  boundary  layer,  and  therefore  shock  wave  boundary  layer  interactions  are 
more  likely  to  affect  upstream  flow  in  a  quiet  tunnel  which  has  a  laminar  boundary  layer  on 
the  nozzle  wall.  If  the  disturbance  is  too  large,  the  tunnel  will  not  start  properly  (f.e.,  will 
not  achieve  uniform  flow  at  the  design  Mach  number  upstream  of  the  test  model). 

The  second  objective  of  the  research  is  to  evaluate  the  capability  of  the  BAM6QT  to  start 
in  the  presence  of  test  models  of  different  bluntness  and  constant  blockage.  The  test  model 
is  a  cone  with  fixed  base  diameter  (and  hence,  fixed  blockage)  and  different  cone  angles  (and 
hence,  different  shock  strengths).  The  results  are  presented  in  Section  3. 


2  Automated  Design  of  Quiet  Hypersonic  Tunnels 

The  concept  of  automated  optimal  design  of  the  test  section  for  a  quiet  hypersonic  tunnel  is 
illustrated  in  Fig.  2.  First,  an  inviscid  contour  is  designed  using  the  Method  of  Characteristics 
(MOC)  with  the  Sivell’s  code  10].  Second,  the  compressible  boundary  layer  is  calculated 
and  the  inviscid  contour  is  corrected  for  this  boundary  layer  growth.  Third,  the  corrected 
nozzle  contour  is  provided  to  a  grid  generator  and  a  mesh  is  created  within  the  test  section. 
Fourth,  a  Navier-Stokes  code  solves  the  steady-state  flow  using  the  grid.  Fifth,  this  flowfield 
is  input  to  a  stability  program  that  computes  the  most  unstable  disturbance  and  establishes 
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the  point  of  transition  along  the  nozzle  contour  using  an  assumed  growth  factor  N.  This 
transition  point  is  used  to  evaluate  the  objective  function  which  is  the  length  of  the  quiet 
zone  as  measured  on  the  test  section  centerline.  This  length  is  defined  by  the  back-to-back 
cones  created  by  the  inviscid  Sivell’s  code  (the  forward-facing  cone)  and  the  Mach  cone 
emanating  from  the  transition  location  (the  backwards- facing  cone)  Lastly,  the  optimizer 
uses  this  value  to  create  a  response  surface  which  is  used  to  select  a  new  set  of  parameters 
to  use  in  the  inviscid  contour  design. 


flow  fold 


final  noxil* 
contour 


GddPro 

Figure  2:  Automated  Optimization  Loop 


The  geometrical  domain  evaluated  by  the  Sivells  code  is  shown  in  Fig.  3.  The  computation 
starts  at  the  throat  using  a  transonic  perturbation  scheme.  The  transonic  perturbation 
scheme  is  valid  near  Mach  1  and  requires  the  throat  radius  of  curvature.  In  this  study, 
the  BAM6QT  coordinates  are  used  for  the  contraction  section  upstream  of  the  throat,  and 
the  throat  radius  of  curvature  specified  in  the  Sivells  input  is  smoothly  connected  to  the 
upstream  contraction  section. 

In  addition  to  the  throat  radius  of  curvature,  the  input  to  the  Sivells  code  includes  the  Mach 
number  along  the  nozzle  centerline  for  some  distance  (until  the  downstream  parallel  flow 
occurs).  The  code  determines  the  appropriate  centerline  Mach  number  distribution  that 
would  result  in  a  nozzle  wall  with  continuous  wall  slope  and  curvature  from  the  throat  to 
the  exit.  The  centerline  Mach  number  distribution  is  based  on  polynomial  functions  that 
match  the  flowfield  solutions  at  each  flow  region  boundary.  Once  this  distribution  has  been 
specified,  a  MOC  procedure  is  used  in  an  inverse  design  mode  to  determine  an  inviscid  nozzle 
wall  contour  that  produces  the  desired  uniform  exit  flow.  This  inviscid  contour  can  then  be 
scaled  to  any  particular  size,  based  on  a  specified  throat  radius  or  exit  radius. 
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Figure  3:  Radial  Flow  and  Expansion  Regions  of  Inviscid  Contour [16] 


The  Sivells  code  is  capable  of  generating  nozzles  with  a  region  of  radial  flow  between  the 
initial  expansion  region  and  the  inflection  point,  where  the  wall  is  shaped  to  turn  the  initial 
characterises  so  as  to  produce  a  uniform  exit  flow.  This  feature  can  be  used  to  suppress  the 
Goertler  instability. 

In  summary,  the  Sivells  code  requires  input  of  four  parameters:  RC  (the  ratio  of  throat 
radius  of  curvature  to  the  throat  radius),  77  (the  expansion  angle  at  the  inflection  points  G 
and  A  in  Fig.  3),  bmach  (the  Mach  number  at  point  D )  and  the  design  test  section  Mach 
number  (taken  to  be  M  =  6  in  this  study). 

The  correction  to  the  contour  due  to  the  compressible  displacement  thickness  of  the  laminar 
boundary  layer  is  calculated  using  using  EDDYBL  19].  A11  example  of  the  correction  is 
shown  in  Fig.  4. 


s 

it  0.25 
1 
■g 

i  0.2 


Inviscid 
Final  Contour 


0.1 


X  Coordinate  (Jt) 


Figure  4:  Example  of  Boundary  Layer  Correction  for  Mach  6  Test  Section 

The  commercial  software  Gridpro  [21  is  used  to  generate  the  structured  grid.  GridPro  is  a 
general  purpose  elliptic  mesh  generator  that  creates  a  three  dimensional,  multi-block  struc- 
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turcd  grid  by  a  variationally-base  method  with  ail  iterative  updating  scheme.  The  user  input 
for  Gridpro  are  the  surface  geometry  and  the  block  topology.  It  is  important  to  note  that 
the  same  topology  definition  is  used  during  the  optimization  while  the  surface  definition 
changes.  An  example  of  the  grid  topology  is  shown  in  Fig.  5 


X  Coordinate  (m) 


Figure  5:  Example  of  Structured  Grid 

The  flow  solver  is  GASPex  [3]  (General  Aerodynamic  Simulation  Program),  a  structured 
multiblock  CFD  solver  for  compressible  viscous  flow.  The  two  dimensional  axisymmetric 
Navier  Stokes  equations  are  solved  using  a  finite  volume  spatial  discretization.  The  opti¬ 
mization  simulations  use  the  steady-state  formulation.  The  computations  are  second  order 
accurate  both  in  time  and  space. 

The  stability  analysis  is  performed  using  the  linear  Parabolized  Stability  Equations  (PSE) 
coupled  with  the  eA  method  to  predict  the  location  of  laminar-to-turbulent  transition.  The 
code  is  LASTRAC  [8],  Since  the  PSE  constitute  an  initial  boundary  value  problem,  the 
wavenumber  and  eigenfunction  at  the  starting  location  are  obtained  from  linear  stability 
theory  (LST).  A  series  of  starting  locations  and  frequencies  are  evaluated,  and  the  growth 
factors  calculated.  An  example  is  shown  in  Fig.  6.  The  furthest  upstream  location  at  which  a 
particular  disturbance  reaches  N  =  6  is  defined  as  the  beginning  of  transition.  The  stability 
analysis  incorporates  the  Goertler  instability  [8].  A  separate  a  posteriori  stability  analysis 
is  performed  to  examine  the  first-  and  second-mode  instabilities  using  the  STABL  code  [7]. 
It  is  observed  that  the  Goertler  instability  dominates  transition  in  nearly  every  case. 

The  computed  length  of  the  quiet  zone  in  the  test  section  (formed  by  the  back-to-back 
cones)  is  the  merit  function  of  the  optimization,  and  the  objective  is  to  determine  the  input 
parameters  RC ,  rj  and  bmach  for  the  Sivells’  code  that  maximize  the  length  of  the  quiet 
zone.  The  test  section  Mach  number  is  fixed  at  M  =  6.  To  improve  the  efficiency  of 
the  optimization,  the  individual  results  for  each  test  section  shape  are  parameterized  as  a 
Response  Surface  [10]  which  is  used  by  the  Multi-Objective  Genetic  Algorithm  (MOGA) 
search  engine  to  determine  the  optimal  design.  The  optimization  is  performed  using  the 
modeFrontier  [1]  environment.  The  workflow  diagram  is  shown  in  Fig.  7. 
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Figure  6:  Example  of  Stability  Analysis 
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Figure  7:  Workflow  Diagram  for  Optimization 


An  example  of  an  optimal  design  for  a  Mach  6  quiet  tunnel  is  shown  in  Fig.  8  with  the 
corresponding  stability  analysis  using  LASTRAC  and  STABL  in  Fig.  9.  A  unique  feature  is 
the  compound  curvature  of  the  nozzle  which  results  in  a  long  section  of  quiet  flow.  This  shape 
is  unconventional,  and  illustrates  the  advantages  of  automated  optimization  (compared  to 
trial-and-error)  to  discover  new  designs.  The  stability  analyis  for  this  design  confirms  that 
the  transition  is  dominated  by  the  Goertler  instability. 
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Figure  8:  Optimal  Design  for  Mach  G  Quiet  Tunnel  (vertical  scale  enlarged) 
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Figure  9:  Stability  Envelope  for  Optimal  Mach  6  Design 

3  Evaluation  of  Tunnel  Startup  for  BAM6QT 

The  BAM6QT  was  modified  [9]  in  2006-2007  to  incorporate  an  expanded  test  section  (Fig.  10) 
for  the  purpose  of  testing  larger  and  blunter  models.  In  the  original  (unmodified)  test  section 
design  (Fig.  1),  the  shock  wave  generated  by  the  test  model  interacts  with  the  laminar 
boundary  layer  on  the  nozzle  wall.  The  concept  of  the  expansion  test  section  is  to  mitigate 
the  shock  wave  laminar  boundary  layer  interaction  on  the  nozzle  wall  by  substituting  a  shock 
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wave  shear  layer  interaction.  The  abrupt  expansion  of  the  test  section  eauses  a  separation  of 
the  incoming  laminar  boundary  layer  and  the  formation  of  a  recireulation  region  and  shear 
layer.  The  principal  question  is  whether  or  not  this  coneept  would  enable  the  tunnel  to  start 
with  larger  test  models. 
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Figure  10:  Expansion  Section  for  Boeing-AFOSR  Mach  6  Quiet  Tunnel 

The  objective  of  the  research  is  to  examine  the  tunnel  startup  process  for  a  test  model  in 
the  expansion  test  section.  A  sharp  cone  at  zero  angle  of  attack  was  selected  for  the  study. 
The  eone  geometry  is  characterized  by  the  base  diameter  and  cone  angle.  The  base  diameter 
determines  the  blockage  ratio  (teM  the  ratio  of  projected  model  cross-seetional  area  to  tunnel 
cross-sectional  area),  while  the  cone  angle  affeets  the  strength  of  the  shock  wave  generated 
by  the  model.  The  bloekage  ratio  lias  traditionally  been  used  to  estimate  the  maximum 
allowable  size  for  a  test  model  To];  however,  it  is  widely  recognized  that  the  results  on 
supersonic  tunnel  starting  do  not  correlate  simply  with  bloekage  due  to  the  differences  m  the 
shoek-boundary  layer  interaction  generated  by  different  models.  Since  there  is  no  universally 
acceptable  criteria  for  determining  the  effeet  of  model  size  on  tunnel  starting  for  supersonie 
and  hypersonie  tunnels,  a  fixed  bloekage  ratio  (r.c.,  fixed  diameter)  was  eliosen  for  simplicity. 

A  7°  half-angle  eone  with  a  5.5  inch  base  diameter  is  the  largest  model  that  ean  be  started 
in  the  BAM6QT  unmodified  test  section  at  zero  angle  of  attaek.  Therefore,  a  fixed  5.5  inch 
base  diameter  was  selected.  The  research  examined  cone  half-angles  ranging  from  15°  to  75° 
to  determine  the  effeet  of  the  shock  strength  on  tunnel  startup. 

A  total  of  six  computations  were  performed.  One  simulation  w*as  performed  for  the  empty  ex¬ 
pansion  test,  seetion,  and  five  simulations  were  performed  for  cone  half-angles  of  15°,  20°,  30°,  50°, 
and  75°.  The  empty  expansion  test  seetion  and  the  15°  half-angle  cone  have  a  10°  compres¬ 
sion  eorner  at  the  downstream  end  of  the  test  seetion,  but  all  other  computations  use  a  more 
gradual  1°  compression  based  upon  subsequent  discussions  with  Prof.  Steve  Sehneider. 

The  computational  grids  were  generated  with  GridPro  [2].  The  domain  for  the  empty  seetion 
consists  of  one  zone  with  23,760  eells  (100  x  241)  and  the  domains  for  the  cone  eases  consist 
of  four  zones  and  contain  approximately  39,000  points.  Grid  clustering  is  performed  with 
a  stretehing  parameter  of  1.105  and  a  first  cell  height  of  10“4  ft  along  the  wall  and  cone 
surfaces  in  order  to  resolve  the  boundary  layers.  Time-aeeurate  simulations  are  performed 
with  GASPex  Version  4.1.0  [3],  using  an  implicit  dual  time  stepping  algorithm  with  a  time 
step  of  10-6  s  and  10  inner  cycles.  Each  case  runs  until  unstart  oceurs  or  until  it  appears 
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steady  (10  —  20  ms  physical  time).  The  inviscid  fluxes  are  discretized  using  Roe’s  scheme  in 
conjunction  with  the  Minmod  limiter  and  MUSCL  spatial  reconstruction.  The  computations 
are  second  order  accurate  both  in  time  and  space.  The  inflow  boundary  condition  is  a  user- 
specified  flow  with  a  boundary  layer  thickness  corresponding  to  a  laminar  boundary  layer 
that  has  been  developing  from  the  stagnation  point  at  the  bleed  lip  near  the  tunnel  throat, 
and  a  uniform  Mach  6.1  freestream  flow.  The  boundary  layer  profile  contains  50  points  for 
the  empty  section  and  40  points  for  the  sections  with  cones,  with  a  boundary  layer  thickness 
of  0.6  in  (15  mm).  The  profile  is  calculated  with  EDDYBL  19],  assuming  a  stagnation 
pressure  of  90  psia,  stagnation  temperature  of  433  K  and  choked  flow  at  the  throat. 

The  initial  conditions  are  chosen  to  approximate  the  actual  tunnel  startup  process.  Warm- 
brod  and  Struck  [18]  describe  the  flow  development  in  a  Ludweig  tube  as  a  two  stage  process 
consisting  of  an  expansion  wave  followed  by  a  shock  wave.  The  rupture  of  the  diaphragm 
(dowustream  of  the  test  section)  generates  an  expansion  fan  that  propagates  upstream.  This 
generates  an  initial  downstream  flow  in  the  test  section.  When  the  expansion  fan  reaches 
the  tunnel  throat,  the  flow  at  the  throat  becomes  sonic  and  an  approximately  normal  shock 
wave  is  formed  in  the  diverging  section  of  the  nozzle  and  propagated  downstream.  The  shock 
eventually  moves  to  the  exit  of  the  tunnel,  accelerating  the  flow  to  Mach  6  within  the  test 
section.  The  initial  conditions,  along  with  the  inflow  conditions,  create  a  shock  wave  that 
propagates  into  the  computational  domain  in  analogy  to  the  actual  tunnel  startup  process. 
The  flow  behind  the  shock  wave  is  the  proper  tunnel  design  flow,  and  is  determined  by  the 
inflow  boundary  condition.  Thus,  while  the  initial  conditions  do  not  precisely  simulate  the 
actual  tunnel  starting  phenomena,  they  qualitatively  model  the  second  stage  of  the  startup 
process  (in  which  a  shock  wave  travels  through  the  test  section  accelerating  the  flow  to  design 
conditions)  and  provide  the  proper  test  section  inflow  conditions  for  the  started  tunnel. 

The  steady  state  (more  precisely,  the  statistically  stationary  state)  flowfield  in  the  empty 
expansion  test  section  is  shown  in  Fig.  11.  The  numerical  Schlieren  is  shown  at  t  =  20  ms. 
The  dark  lines  indicate  a  density  gradient  ie.,  shock  wave,  expansion  fan  or  shear  layer. 
The  expansion  test  section  generates  a  complex  flowfield.  A  conical  expansion  fan  occurs  at 
the  45°  corner  and  intersects  the  centerline.  The  recirculation  region  is  separated  from  the 
core  flow  by  a  shear  layer.  Numerous  weak  shocks  radiate  from  this  shear  layer  and  coalesce 
along  the  centerline  downstream  of  the  expansion  section. 

Figs.  12  and  13  display  the  numerical  Schlieren  image  and  Mach  contours  of  the  flowfield  in 
the  test  section  with  a  15°  half-angle  cone  at  t  =  20  ms  (ie.,  at  the  statistically  stationary 
state).  A  shear  layer  appears  at  the  45°  expansion  corner  and  a  series  of  shocks  turn  the 
flow  back  at  the  10°  compression  corner.  A  recirculation  region  is  formed  in  the  wake  of  the 
cone  and  in  the  expansion  corner.  The  flowfield  upstream  of  the  cone  is  uniform  (outside 
the  boundary  layer).  It  is  therefore  evident  that  the  tunnel  has  successfully  started. 

Contours  of  streamwise  velocity  are  shown  in  Fig.  14  for  the  20°  half-angle  cone  at  t  =  20  ms. 
The  tunnel  appears  to  have  successfully  started  for  this  configuration.  However,  there  is  a 
very  small  upstream  influence  whereby  the  separation  bubble  spreads  beyond  the  expansion 
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Figure  11:  Numerical  Sohlieren  of  Empty  Expansion  Test  Section 
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Figure  12:  Numerical  Selilieren  of  Expansion  Test  Section  with  15°  Cone 


Figure  13:  Mach  Contours  of  Expansion  Test  Section  with  15°  Cone 

corner.  Unlike  the  blunter  cones  (see  below),  this  separation  bubble  is  steady  and  never 
reaches  the  upstream  boundary. 

The  numerical  Selilieren  images  for  cone  half- angles  of  30°,  50°  and  75°  at  3  ms  are  shown  in 
Figs.  15  to  17.  In  each  case  a  separation  bubble  forms  at  the  expansion  corner  and  extends 
upstream.  The  separation  bubble  continues  to  grow  until  it  reaches  the  inflow  boundary. 
The  streamwise  velocity  at  the  inflow  boundary  becomes  negative.  At  this  point,  the  results 
cease  to  be  meaningful  and  the  computation  is  stopped.  These  results  imply  that  the  tunnel 
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Figure  14:  Streamwise  Velocity  Contours  of  Expansion  Test  Section  with  20°  Cone 


cannot  start  with  the  30°,  50°  and  75°  half-angle  cones. 
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Figure  15:  Numerical  Schlieren  of  Expansion  Test  Section  with  30°  Cone  at  t  =  3  ms 


Figure  16:  Numerical  Schlieren  of  Expansion  Test  Section  with  50°  Cone  at  t  =  3  ins 
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Figure  17:  Numerical  Schlieren  of  Expansion  Test  Section  with  75°  Cone  at  t  =  3  ins 
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